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ABSTRACT 

We present a new model atmosphere analysis of the most metal contaminated white dwarf known, 
the DBZ SDSS J073842. 56+183509. 06. Using new high resolution spectroscopic observations taken 
with Keck and Magellan, we determine precise atmospheric parameters and measure abundances of 
14 elements heavier than helium. We also report new Spitzer mid-infrared photometric data that are 
used to better constrain the properties of the debris disk orbiting this star. Our detailed analysis, 
which combines data taken from 7 different observational facilities (GALEX, Gemini, Keck, Magellan, 
MMT, SDSS and Spitzer) clearly demonstrate that J0738+1835 is accreting large amounts of rocky 
terrestrial-like material that has been tidally disrupted into a debris disk. We estimate that the body 
responsible for the photospheric metal contamination was at least as large Ceres, but was much drier, 
with less than 1% of the mass contained in the form of water ice, indicating that it formed interior to 
the snow line around its parent star. We also find a correlation between the abundances (relative to 
Mg and bulk Earth) and the condensation temperature; refractory species are clearly depleted while 
the more volatile elements are possibly enhanced. This could be the signature of a body that formed 
in a lower temperature environment than where Earth formed. Alternatively, we could be witnessing 
the remains of a differentiated body that lost a large part of its outer layers. 
Subject headings: planetary systems - stars: abundances - stars: atmospheres — white dwarfs 



L INTRODUCTION 

The heavy elements observed in the spectra of other- 
wise pure hydrogen or helium atmosphere white dwarfs 
represent handy fingerprints of asteroids/rocky plan- 
ets that survived the late phases of s tellar evolu- 
tion (iZuckerman et al.l [20071 : iKlein et al.l IMl I20TTI : 
IDufour et al.ll2010n . There is growing evidence that per- 
haps all metal polluted white dwarfs of spectral type 
DZ, DBZ and DAZ have acquired their heavy mate- 
rial from an orbiting debris disk reservoir whose ori- 
gin is explained by the tidal disruption of one or many 
large rocky bodies that ventured too close to the star 
jDebcs & Sigurdsson 2002; Jura 20 0"1 l2006t lJura et al.l 
[2007l:lJuraii2008l:iFarihi et al.li2010aiibr[Melis et al.ll2010D . 
Such disks, which are easily detectable at infrared wave- 
lengths, are now discovered with an accelerating pace 
with more than 20 cases unc o vered in the l ast 6 years 
alone (see iFarihi et al.l 120091 : iFarihil 120111 : iKilic et al.l 
I2OIII and references therein). 

Meanwhile, high resolution observations in the optical 
also revealed the presence of weak calcium lines in many 
white dwarfs, suggesting that at least 25% of their pro- 
genitors possessed orbiting planetary bodies with masses 
often co mparable to that of the larges t solar system as- 
teroids {ZuctermaEiEin [IQOI [20l3) . Unfortunately, 
the majority of the known metal polluted white dwarfs 
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only show one or two heavy elements in the optical, 
hardly en ough to make deta iled planetary composition 
inquiries (jDufour et al.l l200"7[ ) . In a few extreme cases, 
however, a handful of elements are observed, opening 
a unique window for studying rocky exoplanet compo- 
sitions. The first investig ations, those of GD 362 and 
CD 40 (jZuckerman et al.i [20071 IKlein et al.l [20Toh re- 
vealed striking similarity between the accreted bodies 
and bulk Earth material, providing the first comprehen- 
sive measurements of the bulk elemental composition of 
ancient planetary systems. Analysis of similar objects, 
based on Keck high resolution observations, quickly fol- 
lowed these pionering studies, bringing to 6 the num- 
ber of metal polluted white dwarfs with more than 8 
heavy elements detected in the optical (the previous two 
along with G2 41-6, NLTT 43806, PG 1225-079 and 
HS 2253-l-8023. [Zuckerman etliI1[2010l [20Tlt IKlein et al.l 
[2011). 

Although the global abundance patterns found for 
these few objects were relatively well explained by the ac- 
cretion of rocky planetary material similar to bulk Earth, 
closer looks reveal a wide range of compositional diversity 
from star to star (or, equivalently, planets to planets), 
especially for the most refractory elements . Interest- 
ingly, d ynamical simulations such as those of iBond et al.l 
(|2010al |bl) also predict a wide variety of chemical prop- 
erties and structure compared to what is observed in 
the solar system. The most metal polluted white dwarfs 
could thus represent perfect testbeds to discriminate be- 
tween various planet formation scenarios, allowing the 
first empirical verification of extrasolar terrestrial planet 
formation simulations. However, other processes such 
as post AGB thermal heating, collisions and wind strip- 
ping, could also play an important role in the final bulk 
composition that is measured at a white dwarf photo- 
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sphere. Understanding the physical mechanisms respon- 
sible for the various observed elemental abundance pat- 
terns represents one of the biggest challenge of this rela- 
tively young field of research. As the sample of well stud- 
ied metal contaminated white dwarfs increases, a more 
comprehensive picture should emerge, leading to a better 
understanding of how, both individually and statistically, 
these plane tary systems form a nd evolve. 

Recently, iDufour et al.l ()201Cl[ ) discovered an extremely 
polluted DBZ white dwarf, SDSS J073842. 56+183509.06 
(hereafter J0738+1835) in the Sloan Digital Sky Sur- 
vey. A preliminary analysis, based on medium res- 
olution optical spectroscopy, revealed the presence of 
O, Mg, Si, Ca and Fe in record quantities. Given 
that Keck/HIRES spectroscopy of the heavily polluted 
white dwarf GD 362 revealed 15 different heavy ele- 
ments (jZuckerman et al.l l2007f l . three of which are de- 
tecte d in low resolution spectroscopy (jGianninas et al.l 
120041 ) ■ we anticipated that higher resolution observa- 
tions of J0738-I-1835 would reveal several a dditional ele- 
ments not detected in IDufour et al.l (|2010D 's low resolu- 
tion data. 

This paper reports a detailed analysis of J0738+1835 
based on new high resolution optical spectra taken at 
the Keck and Magellan telescopes. We also report new 
mid-infrared photometric data taken with Spitzer which, 
when combined with Gemini JHK photometry, provides 
better constraints on the debris disk physical properties. 
In addition to the heavy elements already detected by 
IDufour et al.l (j2010D , our new spectroscopic observations 
of J0738-I-1835 reveal the presence of 9 new elements 
(Na, Al, Sc, Ti, V, Cr, Mn, Co and Ni), for a total 
of 14 chemical species heavier than helium. This is the 
second highest number of heavy elements ever observed 
at a white dwarf photosphere ( the champion, w'ith 15 
different metals, being GD 362: iZuckerman et al.l[2007h 
thus making J0738-I-1835 a new and important member 
of the very select group of polluted white dwarfs with 
eight or more elements heavier than helium detected op- 
tically (Zuckcrman ct al. 2011). 

In §[21 we describe the new observations. Our detailed 
analysis of the spectroscopic and photometric observa- 
tions follows in § [3] while the results are discussed and 
summarized in § |3] and § [5] 

2. OBSERVATIONS 

Our detailed study of SDSS J0738+1835 uses data 
taken from 7 different observational facilities (GALEX, 
SDSS, Gemini, MMT, Keck, Magellan and Spitzer), with 
a strong emphasis on new Keck HIRES spectroscopy. 
Part of these data ha ve already been pres ented in the 
preliminary analysis of IDufour et al.l ()2010f ) and are not 
discussed further here. Below we describe the new data. 

2.1. Keck/HIRES 

We us ed the Hig h Resolution Echelle Spectrometer 
(HIRES, iVogt et al.l[l994l ) on the Keck I 10m telescope 
(Mauna Kea Observatory) to obtain high resolution op- 
tical spectroscopy of J0738-I-1835 over 4 hours on UT 
2011 March 23. We use the blue collimator and the 
1.148" slit with the cross disperser angle of 1.075, pro- 
viding wavelength coverage from 3100 A to 5990 A with 
a resolving power of 37,000. The seeing started at 0.9" 
with an average around 1.1 ". We use the MAuna Kea 



Echelle Extraction (MAKEE) software for data analysis. 
We flat-field the exposures, optimally extract the spec- 
tra from the two dimensional image traced by the stan- 
dard star Feige 34 observed on the same night, remove 
cosmic rays, and wavelength calibrate the spectra using 
Th-Ar lamps. MAKEE does not extract spectra in a few 
of the partially covered orders. We use the HIRES Re- 
dux package written by J. X. Prochaska to extract these 
missing orders from the MAKEE pipeline. The HIRES 
observations of J0738-I-1835 is thus composed of 58 short 
spectroscopic segments, typically 60 to 90 A wide, with 
some wavelength overlap in adjacent orders. The over- 
lapping parts have been co-added to slightly improve the 
signal-to-noise ratio. 

2.2. Magellan/MagE Optical Spectroscopy 

Moderate-resolution optical spectra of J0738+1835 
were obtained on 201 1 March 19 (UT) wi th the Magel- 
lan Echellette (MagE: lMarshall et al]l2008f ). mounted on 
the 6.5m Landon Clay Telescope at Las Campanas Ob- 
servatory. Conditions during the observations were clear 
with 0.6" seeing. Two exposures of 1200 s each were ob- 
tained over an airmass range of 1.48-1.49 using the 0.7" 
slit aligned with the parallactic angle; this setup pro- 
vided 3200-10050 A spectroscopy at a resolution « 8000 
(although only the red part longward of 6000 A, where 
some important oxygen lines are present, is used in the 
present analysis since the Keck observations covers the 
rest with a better resolution and signal-to- noise ratio). 

We also obser ved the spe ctrophotometric flux stan- 
dard Hiltner 600 (jHamuv et al. 1994) on the same night 
for flux calibration. Th-Ar lamps were obtained after 
each source observation for wavelength calibration. Or- 
der tracing and pixel response were calibrated with inter- 
nal Xe lamp flats and twilight sky flats, respectively, ob- 
tained at the beginning of the night. D ata were reduced 
using the MASE reduction pipeline ([Bochanski et al.l 
|2009( ). following standard procedures for order tracing, 
flat field correction, vacuum wavelength calibration (in- 
cluding heliocentric correction), optimal source extrac- 
tion, order stitching, and flux calibration. 

2.3. MMT spectroscopy 

We used the 6.5m MMT equipped with the Blue 
Channel spectrograph to obtain medium resolution spec- 
troscopy of J0738+1835 over 109 minutes on UT 2010 
March 19. We operate the spectrograph with the 832 
line mm^^ grating in second order, providing wavelength 
coverage from 3170 A to 4100 A and a spectral resolution 
of « 1 A. All spectra were obtained at the parallactic an- 
gle. We used He-Ne-Ar comparison lamp exposures and 
blue spectrophotometric standards (|Massev ei~aII[T988[ ) 
for wavelength and flux calibration, respectively. We re- 
duce the data using standard IRAF routines. 

2.4. Spitzer 

We used the warm Spitzer equipped with the InfraRed 
Array Camera (IRAC, Fazio et al. 2004) to obtain mid- 
infrared photometry of J0738-hl835 on UT 2010 Decem- 
ber 1 as part of program 70023. We obtained 3.6 and 
4.5 ^m images with integration times of 100 seconds for 
nine dither positions. We use the IDL astrolib pack- 
ages to perform aperture photometry on the individual 
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corrected Basic Calibrated Data fraraes from the latest 
available pipeline reduction. We get similar results using 
2 and 3 pixel apertures, but we quote the results using 
the smallest aperture since it has the smallest errors. 

Following the IRAC calibration procedures, we cor- 
rect for the location of the source in the array before 
averaging the fluxes of each of the dithered frames at 
each wavelength. We also correct the Channel 1 (3.6 
^m) photometry for the pixel-phase-dependence. We 
estimate the photometric error bars from the observed 
scatter in the nine images corresponding to the dither 
positions. We add the 3% abso lu te calibratio n error 
in quadrature ([Reach et al.ll2005[ ). iReach et all (|2OO90 
demonstrate that the color corrections for dusty white 
dwarfs like G29-38 are smaU (0.4-0.5%) for channels 1 
and 2. We ignore these corrections for J0738-t-1835. 

The final observed flux and uncertainties at 3.6 fim and 
4.5 fim are 81.1 ± 2.7 iiJy and 76.7 ± 2.7 fiJy respec- 
tively. 

3. DETAILED ANALYSIS 
3.1. Model Atmospheres and Fitting Technique 

The model atmosphere and synthetic spectrum calcu- 
lations are performed using the same code that was used 
in th e flrst attempt to model J0738-I-1835 (jPufour et al.l 
l2010f l with the exception that we now use atomic data 
provided by the Vienna Atomic Line Database (VALEEI) 
instead of the Kurucz line list. 

As a first step, we calculated a model structure us- 
ing the same param eters (Toff, logg and abundances) as 
iDufour et al.l (|2010f ). including all the transitions (from 
200 to 10,000 A) present in our linelist. Using this ther- 
modynamic structure, we next calculated several grids of 
synthetic spectra, one for each element of interest, keep- 
ing the abundance of all other elements fixed to their 
original values. For example, the Fe grid covers a range 
of abundance from log [n(Fe)/n(He)] — —4.0 to —6.0 in 
steps of 0.5 dex with the abundances of all other met- 
als fixed to those previously determined. The various 
grids typically vary by log [n(Z)/n(He)] — ±1.5 dex (in 
step of 0.5 dex) around the previously deter mined (or 
expec ted by assuming CI chondrites ratios from lLoddersI 
120031) abundance. These grids are then used to determine 
the different elemental abundances as follows. 

Since the optical spectra for J0738-I-1835 are extremely 
crowded with several hundred lines, many of which over- 
lap, it is not practical to fit each line individually using 
equivalent width measurements. Instead, we made short 
segments, typically 15 to 30 Angstrom wide, centered on 
each of the strongest absorption features observed in our 
dataset (i.e., those that are unambiguously detected and 
deeper than ~10% of the continuum flux level). Then, 
each segment (a few hundreds) arc fltted separately with 
the appropriate grids. 

Our fits are done by minimizing the value of the 
taken as the sum of the difference between the observed 
and model fiuxes over the wavelength range of interest, 
all frequency points being given an equal weight. We 
consider the abundance, the solid angle and a linear term 
to account for the slope (the numerous HIRES orders 
are not fiux calibrated) as free parameters. This method 

^ http:/ /vald. astro. univic.ac.at/~vald/php/vald.php 



allows a very good determination of the local continuum 
and naturally takes into account the contribution of the 
absorption from other elements. It is thus particularly 
well suited for situations where it is difficult to precisely 
determine equivalent widths of individual lines, either 
because line blending is important or the signal-to-noise 
ratio is marginal. We show in Figure 1 typical examples 
of segment fits. In order to better visualize the fitted 
features, we plotted in blue a synthetic spectrum in which 
the element that is fitted has completely been removed. 
We note that in many cases, more than one line from a 
given fitted element can be found in a single segment. In 
these cases, the determined abundance of the segment is 
essentially a weighted average of the lines, more weight 
(more frequency points) being put on the strongest lines. 
We finally take as the measured abundance of an element 
the average over all the segments. 

It must be noted, however, that these abundances 
were determined from synthetic sp ectra calculated with 
a structure corresponding to the IDufour et al.l ()2010f ) 
solution. The presence of new species, as well as the 
slightly different abundances that we now derive (and, to 
a smaller extent, the use of the VALD data), has a small 
impact on the thermodynamic structure of our model. 
Hence, in order to obtain a self-consistent solution, it 
is necessary to repeat the above mentioned procedure 
(i.e., recalculate all the grids and redo the fits) but, this 
time, starting from a model structure calculated with 
the newly determined set of abundances. We iterate un- 
til the input abundances in the model structure and the 
final measurements from our grids are in agreement. 

However, even this solution can not be considered fi- 
nal since such an analysis is done for a fixed value of Teg 
and log 5 (13,600 K and 8.5, respectively). These pa- 
rameters were originally determined from fitting the Hel 
line profiles as well as simultaneo usly accommodatin g the 
Mgl and M gll line streiigths (seelDufour et al.ll2010D . As 
discussed in lKlein et al. I (|2oIa 120117 the absolute abun- 
dance of the various elements can vary significantly (up 
to 0.4 dex in some stars) with small variations of Toff and 
log 5. The relative composition of the polluting elements, 
however, are much less sensitive to the exact final param- 
eters adopted and can thus be used as a high precision 
measurement of the accreted planets/planetesimals (see 
below). Nevertheless, we believe it is best to take full 
advantage of the opportunities provided by our dataset, 
and proceed in obtaining the best atmospheric parame- 
ters as possible by reevaluating the effective temperature, 
the surface gravity and the elemental abundances in a self 
consistent way. 

We thus repeated the whole analysis as described above 
but, this time, starting with model structures calculated 
at effective temperatures of 13,300 and 13,900 K (cor- 
responding to the estimated uncertainties of ±300 K in 
IDufour et al.ll2010l ). We then focus on elements for which 
lines from two different ionization state are present. In 
particular, we select a subsample of the strongest well 
isolated lines of Mg and Fe, the elements that show the 
largest number of lines from both the neutral and ion- 
ized state. In Figures 2 and 3, we plot the average abun- 
dance for each ion as a function of effective tempera- 
ture, taking the standard deviation of the set of indi- 
vidual abundances as error bars. Clearly, a much better 
agreement between the two ionization states is obtained 
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Fig. 1. — Typical example of fits to small HIRES segments. Only one element at a time is fitted in eacfi panel, all other elements being 
kept fixed to the value determined in the previous iteration (see text). The determined abundance of the element of interest, which is not 
the final abundance reported in Table 2, is indicated in the lower left corner of each panel. The red line is our best fit while the blue line 
is a synthetic spectrum where the element of interest has been removed to help localize the regions where the fitted lines contribute. The 
observations have been artificially shifted to air wavelengths in order to facilitate line identification using line lists. 

for a slig htly higher eff e ctive temperature than first ob- 
tained bv lDufour et al.l (|201Cl f). A remarkable agreement 
between the various ion abundances is reached for an 
effective temperature of ^ 13,950 K. This effective tem- 
perature estimate, which is determined solely from Mg 
and Fe lines, is completely independent of the uncertain 
treatment of van der Waals Hel line broadening. 

Next, we take this newly determined effective temper- 
ature and repeat again the above procedure but this time 
with surface gravities between log (7= 8.3 and 8.7. Our 
results, summarized in Figure 4 and 5, indicate that a 
surface gravity lower than log 17= 8.5 better reproduces 
the different ion abundances. A surface gravity of log (7~ 
8.4 is preferred by the Mg lines while a slightly lower sur- 
face gravity is required for the Fe lines, although a value 
of log g^SA is compatible with the error bars. 

We thus adopt, as our final parameters, an effective 
temperature of 13,950 K and a surface gravity of logg« 
8.4 with conservative uncertainties of ±100 K and ±0.2 
dex respectivelj0. These uncertainties are next prop- 
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® Such a level of accuracy for the effective temperature (~0.7%) 
is among the most precise ever determined for a white dwarf star. 
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Fig. 2. — Average abundances of magnesium lines as a function of 
effective temperature for two different ionization states. The error 
bars represent the standard deviation of the set of individual. A 
much better internal coherence is found for an effective temperature 
near 13,950 K. 

erly propagated for the determination of errors on the 
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Fig. 3. — Same as Figure 2, but for iron. 
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Fig. 4. — Average abundances of magnesium lines as a function 
of surface gravity for two different ionization states. The error bars 
represent the standard deviation of the set of individual. A shghtly 
lower surface gravity than log 3= 8.5 appears to be preferable. 
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Fig. 5. — Same as Figure 4. but for iron. 



cooling age, luminosity, mass, radius and distance of the 
star (obtained using the solid angle from ugriz photom- 
etry). The evolutionary models used are similar to those 



TABLE 1 

Stellar parameters adopted for 
SDSS .70738+1835 



Parameter 


Value 


Teff(K) 


13950 ± 100 


log 9 


8.4 ± 0.2 


Mwd/Mq 


0.841 ± 0.131 


Minit/M0 


4.47 ±0.36 " 


R/Rq 


0.00958 ± 0.0015 


log L/Lq 


-2.50 ±0.12 


D 


147 pc ±23 


Cooling Age 


477 Myr ± 160 


log(A/Hc/Af*) 


-6.41 +0.56/-0.28 



^ Initial mass of the main sequence pro- 
genitor calculated using the Initial-Final 
Mass Relation of lWilliams et al.l ll2009f) 



described in iFontaine et al.l ()200lD but with C/0 cores, 
and thickness of the helium and hydrogen layers of re- 
spectively q(He) = 10^^ and q(H) = 10"^'^, which are 
representative of helium-rich atmosphere white dwarfs. 
To evaluate the abundance uncertain ties, we proceed in 
a similar manner to that described in lKlein et al.l (|201Cll ) 
and vary the model Tcff and log (7 by their uncertainties, 
one at a time, determine the average change in abun- 
dances for each parameter, and add them in quadrature 
to the standard deviation of the set of individual abun- 
dances. Our final parameters are presented in Table 1 
and 2 while our best fit models are plotted over the spec- 
troscopic data for the most interesting HIRES and MagE 
orders in Figures 6 to 12. 

We also, a posteriori, verified explicitly that this new 
solution was still compatible with the lower re solution 
MMT spectroscopic data used in iDufo ur et al.l ([2010). 
Figure 13 shows that the slight change in atmospheric pa- 
rameters that we derive here did not have any noticeable 
impact on the quality of the fit to the various Hel lines 
in our "old" MMT data. In fact, we even observe that 
most of the dis crepancies in the p redicted line strengths 
of Fe lines (see IDufour et al]l2010 '. Figure 3) are now al- 
most completely removed. We attribute this result on 
our improvement in the atmospheric parameters as well 
as the use of better atomic line data (VALD). 

It is interesting to note that while the Hel line profiles 
seem well reproduced by our models when observed at 
relatively low resolution, such is not the case at higher 
resolution. Indeed, detailed inspection of strong Hel line 
at AA 3889, 4471, 4713, 4922 and 5876 indicates that 
the synthetic models poorly reproduce the HIRES ob- 
servations (the MagE observations of Hel AA 5876 and 
6678, taken at lower resolution, are more satisfactory). 
Similar problems in Hel line profiles have also been ob- 
served in other cool DBZ white dwarfs (see, in particu- 
lar, the narrow absorption core compo nent in Hel A5876 
for GD 362, GD 16 and PG 1225-079. IZuckerman et"all 
120071: iKlein et al.l 1201 ID . Note that varying the atmo- 
spheric parameters (Tcff /log (?) from the optimal values 
determined above does not improve our fit to the helium 
lines; the discrepancy appear to reside in the predicted 
line shape. The explanation for the features morphology 
eludes us at this time. Until this issue is resolved, we 
believe it is safer to use the ionization balance of heavy 
elements to obtain the effective temperature and surface 
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TABLE 2 



Element 


log [n{Z)/n{He)]phot 


A/cvz/(102ig) 


log Tset (yr) 


[n{Z)/n(Fe)]acc 


A//(10**g s^i) 


1 H 


-5.73 




0.17 


0.310 


oo 






8 O 


-3.81 




0.19 


407.86 


5.244 


9.52 


740.2 


11 Na 


-6.36 


± 


0.16 


1.639 


5.238 


2.7 X 10-2 


3.02 


12 Mg 


-•4i.Ot5 


m 


U.U 1 


83.33 


5.258 


1.24 


146.38 


13 Al 


-6.39 


± 


0.11 


1.792 


5.244 


2.5 X 10-2 


3.25 


14 Si 


-4.90 


1 

± 


0.16 


57.99 


5.248 


0.77 


104.36 


20 Ca 


-6.23 


± 


0.15 


3.907 


5.044 


5.8 X 10-2 


11.24 


21 Sc 


-9.55 


± 


0.18 


2.05x10"^ 


5.010 


2.9 X 10-5 


6.38x10-3 


22 Ti 


-7.95 


± 


0.11 


8.87x10-2 


5.007 


1.2 X 10-3 


0.278 


23 V 


-8.50 


± 


0.17 


2.65x10-2 


5.006 


3.4 X 10-"' 


8.31x10-2 


24 Cr 


-6.76 


± 


0.12 


1.492 


5.026 


1.8 X 10-2 


4.48 


25 Mn 


-7.11 


± 


0.11 


0.693 


5.028 


7.7 X 10-3 


2.07 


26 Fe 


-4.98 


± 


0.09 


94.91 


5.047 


1.00 


271.32 


27 Co 


-7.76 


± 


0.19 


0.165 


5.042 


1.7 X 10-3 


0.479 


28 Ni 


-6.31 


± 


0.10 


4.721 


5.063 


4.6 X 10-2 


12.997 


Total 
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Fig. 6. — Display of our final solution (in red) over the most interesting HIRES orders. Ticks and labels identify the strongest features 
present in the spectra. The observations have been artificially shifted to air wavelengths in order to facilitate Hne identification using 
available line lists. 



gravity with high precision 

We also present in Figure 14 a similar exercise using 
unpublished wl A resolution MMT spectroscopic obser- 
vations that covers the bluer part of the electromagnetic 



spectrum. While the signal-to-noise ratio and resolution 
of these observations were first judged to be ins ufficient 
to pro vide significant improvements over Dufou ~et al.l 
(|2010D 's solution, they are nevertheless useful to confirm 
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Fig. 7. — Same as previous Figure. 

that no significant change in the elemental abundances 
of the major species has occurred between the ~1 year 
that separates the MMT and Keck observing runs. This 
is not surprising given that this object is most probably 
accreting material at a steady state and that the diffu- 
sion timescales of the heavy elements are of the order of 
10^ years (see below). 

Although the final abundance determination has been 
obtained only from an average of a subsample of all the 
lines present in the Keck and Magellan spectroscopic 
datasets, we find that the synthetic spectrum calculated 
with our final solution agrees remarkably well with the 
observations for the vast majority of the lines over the 
whole spectrum (see Figures 6 to 12). This can most 
probably be attributed to the fact that the method de- 
scribed above has yielded a very accurate set of atmo- 
spheric parameters that minimized the dispersion among 
the set of individual measurements. Indeed, we find that 
the standard deviation of the set of individual measure- 
ments for a given element vary typically between 0.04 to 
0.10 dex (with a few exceptions, notably oxygen, that 
have slightly larger dispersions; see below). 

There are, however, a few individual lines which signifi- 
cantly depart from our final optimal solution (most prob- 
ably due to discrepancies in atomic line data or broad- 
ening parameters). These can be observed in Figures 6 



to 12 as features that are not very well reproduced by 
the final synthetic spectrum (see, for example, lines near 
AA4370, 4740, 4850, 5400, 5505). Overah, fortunately, 
the final atmospheric parameters are not significantly af- 
fected by these lines since their contribution tend, for 
the few included in the final set, to be averaged out 
from using multiple lines from each elements (O, Na, Al, 
Sc and Co, which rely on 3 or less individual measure- 
ments, must thus be considered more uncertain on that 
account). 

We note, finally, that our observations for J0738+1835 
show no sign of emission cores i n the Call H&K 
lines, as ob s erved for PC 1225-079 ijKlein et al.|[2011.) . 
iKlein et al.l (|201lD proposed that NLTE effects could 
possibly cause such emission since these line cores are 
formed high up in the atmosphere. However, comparing 
the thermodyn amic st r uctur es of PG 1225—079, calcu- 
lated with Klei n et al.l (|201l!)'s parameters, with that of 
J0738-fl835 argues against this interpretation (not to 
mention that NLTE effect are not expected in such a 
low effective temperature white dwarf). Indeed, we find 
that the atmospheric pressure in J0738-I-1835 is lower 
than that of PG 1225—079, while its gas temperature is 
higher at every depth point in the atmosphere. Since 
both of these should favor NLTE effects (less collisions 
and a more intense radiation field), if NLTE was the 



8 



X 

N 

• I — I 

a 

o 




0.0 



0.0 



3820 3830 3840 ~ 3850 

Wavelength (A 



3860 



3870 



3880 



Fig. 8. — Same as previous Figure. 

explanation, J0738+1835 (and other similar DBZ white 
dwarfs) would be even more likely to show cores in emis- 
sion. Since they are not observed in any other object, 
PG 1225—079 is more likely a unique and interesting case 
of white dwarf chromospheric activity. 

It is interesting to note, however, that our MagE spec- 
troscopic observations clearly show evidence of the IR 
Call double-peaked emissions (see Figure 12). These 
features, which wer e already observ ed in the SDSS spec- 
troscopic data fsee 'Gansickel [20TTI ). are the typical sig- 
nature of g aseous debris di sks rotating around a white 
dwarf (jG&nsicke'et al.ll2nnd [2?)07l 120081 ) . 

It has been proposed that gas disk-hosting systems 
could be the result of multiple infalling bodies orig inating 
from a solar system-like asteroid belt ()Jurall2008f ). How- 
ever, the total mass of heavy element presently found in 
J0738-|-1835's convective zone represent a very large frac- 
tion of the mass of an asteroid belt. While it is possible 
to explain this large amount by the accretion of several 
small asteroids, as newly-shredded asteroids dust arrives 
near the pre-existing disk, destructive grain-grain col- 
lisions should rapidely cause the disk to evolve toward 
a large gaseous system (Jura 2008). Given that a sig- 
nificant dusty disk component is revealed from a large 
infrared excess (see below), we are most likely withness- 
ing a case where one large object generated the dense 



dust disk which dominate the system and that the small 
amount of gas observed is the result of much smaller 
bodies hitting the dust disk from time to time. It is also 
possible that we are watching the system in a transition 
phase where the dusty disk is starting to slowly dissipate, 
through internal motions, into a gaseous phase. Alter- 
natively, the emission li ne flu x could be understood, as 
proposed by iGansickj ()2011[ ). by some heating of the 
top layers of the disk with ultraviolet p hotons from the 
white dwarf (see also lMelis et al.lf2010( ). Unfortunately, 
our understanding of white dwarfs debris/gaseous disk is 
insufficient at this point to speculate furthermore on the 
formation/evolution of these components. 

3.2. Infrared Excess and Disk Model 

The presence of an accretion disk around J0738-I-1835 
was unamb iguously estab l ished from JHK excess pho- 
tometry in iDufour et al.l ()201ClD . However, due to the 
lack of mid-infrared data, the outer temperature was not 
well constrained. Here we use Spitzer (IRAC) 3.6 and 
4.5 fXTa data combined with Gemini JHK photometry to 
get better constraints on the physical parameters of the 
disk. 

In the same manner as in IDufour et al. 



use the optically thick flat-disk models of 
and determine the best fit disk parameters 



Jural (|2003[ ) 



We obtain 
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Fig. 9. — Same as previous Figure. 

r,„ = 1600 ± 100 K, Tout = QOOljoo K and an inclination 
angle of 40 ± 20 degree (see Figure 15). The uncertainties 
for these par ameters were obta- ined following the method 
described in IKilic et all (j201lD . Briefly, a Monte Carlo 
analysis is performed where the observed photometric 
fluxes f are replaced with f + gSf, where 6i is the error 
in flux and g is a Gaussian deviate with zero mean and 
unit variance. For each of 10,000 sets of modified pho- 
tometry, the analysis is repeated to derive new best-fit 
parameters for the disk. The interquartile range of these 
parameters is adopted as the uncertainty. The larger low 
range uncertainty found for the outer temperature is the 
result of not having data beyond 5 microns since it is 
possible that the disk may extend to larger radii with 
a corresponding cooler lower-limit for the temperature. 
With Tout — 900 K, corresponding to an outer radii of 
'--^ 0.1 R0, the disk is well placed inside the tidal radius 
of the white dwarf (r-^ 1 R©), further strengthening the 
idea that the origin of the disk, as well as the material 
polluting the white dwarf's convection zone, is the tidal 
disruption of a large rocky body. 

4. RESULTS AND DISCUSSION 

As discussed in details in iKoesteii ()2009[ ). the inter- 
pretation of the abundance pattern measured at a white 
dwarf photosphere is not always straightforward. Indeed, 



depending in which of the three phases of the accretion- 
diffusion process the star is observed (which we gener- 
ally do not know, especially in the case of helium-rich 
atmosphere white dwarfs with relatively long diffusion 
timescales) , the relation between the photospheric abun- 
dance and that of the accreted body differs. 

For instance, in the first phase, when the accretion 
process just began, the observed relative abundance of 
each element measured at the photosphere directly trans- 
lates into that of the polluting body. Next, the so-called 
steady state is reached and the abundance pattern can 
no longer be directly interpreted since each element has 
its own characteristic diffusion timescale. Fortunately, 
one can still recover the relative abundance of the source 
material by simply applying a correction that takes into 
account the difference in diffusfo n time of the various 
elements (see lDupuis et allll9 93f). Finally, when the ac- 
cretion ends, each element disappears exponentially with 
their corresponding diffusion constant, quickly affecting 
the abundance ratios between elements with different 
timescales. In the case of J0738-I-1835, given the strength 
of the infrared excess, as well as the large accretion rates 
needed to explain the photospheric pollution (see Table 
2), one can be confident that this system is currently 
accreting and not in the declining phase. 

It is generally assumed that systems with important 
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Fig. 10. — Same as previous Figure. 

infrared excess are in the steady state regime, although 
there is always a small chance that we are observing this 
object in the early phase. The probability of catching 
an object at that moment depends on the disk lifetime, 
which is, ho wever , highl y uncertain. Based on simple 
assumptions, lJural (j2008| ) estimates that a median-sized 
asteroid may produce a disk that last around 1.5 x 10^ 
yrs, comparable to the diffusion timescales of the heavy 
elements in J0738-|-1835's photosphere (see below). If 
dusty disk lifetimes are truly that short, then the steady 
state might not be reached and the observed abundance 
pattern could be directly interpreted. However, it is 
also possible that dusty disks, depending on their vis- 
cosity and composition, may last as long as 10^ yrs, like 
the rings of Saturn. Accordingly, given our poor under- 
standing of the temporal evolution of circumstellar disks, 
caution is well advised when interpreting observed abun- 
dance ratios (for example, see iRafikovi i2011i : ?). The 
presence of a gaseous component, however, seem to indi- 
cate that J0738-|-1835's disk will probably evolve much 
faster than that. Nevertheless, assuming that a much 
larger body such as the one that was tidally disrupted 
near JG738-I-1835 could last longer (a few settling time) 
than the simple estimates of lJural ()2008[ ). then a steady 
state could be reached. 
In any case, the maximum correction applicable on 



the abundance ratio of two elements in the steady state 
regime is only ~1.6 (see below) for J0738-|-1835's atmo- 
spheric parameters, and our conclusions are not signifi- 
cantly affected by the exact regime in which J0738+1835 
is assumed to be. For simplicity, we will thus simply as- 
sume that J0738-f 183 is accreting material in the steady 
state regime for the rest of our discussion. 

4.1. Convection Zone Models and Diffusion Timescales 

The convection zone models for our present analysis of 
J0738+1835 have been built as described in Dufour et al. 
(2010). In brief, we computed full stellar structures spec- 
ified by the values of their surface gravity and effective 
temperature (to ease comparison with the spectroscopic 
observations) , as well as a chemical stratification given by 
a He-rich {Z = 0.001) outer layer of fractional mass log 
(1 — M{r)/M^:) = —3.0 surrounding a C/0 core. Given 
the relative insensitivity of the convection zone thickness 
on the assumed convective efficiency for the atmospheric 
parameters appropriate for J0738-I-1835 (see Dufour et 
al. 2010), we only considered the so-called ML2/q!=1.0 
version in the present paper. Convection zone models 
were thus computed for five values of the effective tem- 
perature (Tefi= 13,000 K to 14,200 K, in steps of 300 
K), and five values of the surface gravity (log 5= 8.1 
to 8.9, in steps of 0.2 dex). This domain in parame- 
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ter space is centered on the atmospheric parameters for 
J0738+1835 initially inferred in Dufour et al. (2010), 
and more than cover the uncertainty range ultimately 
found in the present study. 

The spline coefhcients provided by iPaquette et al.l 
(|l986a) to characterize the collision integrals involved in 
diffusion processes have been widely used in stellar as- 
trophysics. I n the case of white dw arfs, they have first 
been used bv IPaquette et al.l ()1986b[ ) to estimate the dif- 
fusion timescales of a few representative metals at the 
base of the convection zone of sufficiently cool He- and 
H-rich stars. The exact same approach has been followed 
recently by Koester (2009) for additional elements, in- 
cluding the use of the same rough model of pressure ion- 
ization devised bv IPaquette et al.l (|1986b[ ) for estimating 
the average charge of trace heavy elements under white 
dwarf conditions. Given the importance of investigating 
the presence of trace amounts of metals in cool white 
dwarfs in the context of the new paradigm of a plane- 
tary origin, our group at Montreal has gone back to the 
basic problem of computing still more realistic diffusion 
coefficients. The details of this effort will be reported 
elsewhere, but we point out here that 1) improvements 
in the numerical evaluations of the collision integrals have 
been implemented, 2) a more physical description of the 
screening length has been used, 3) an improved model of 



pressure ionization (insuring smoothness of the charge of 
an element with depth) has been adopted, and 4) all 27 
elements from Li to Cu in the periodic table have been 
considered in the calculations. 

Figure 16 summarizes our computations of settling 
timescales at the base of the convection zone for our He- 
rich white dwarf models. In the range of effective tem- 
peratures depicted, the diffusion timescale for a given el- 
ement increases monotonically with decreasing tempera- 
ture at fixed gravity, but the effect is relatively small. On 
the other hand, in the range of surface gravity considered 
in Figure 16, the diffusion timescale increases relatively 
quickly with decreasing surface gravity at fixed temper- 
ature. We note that, while the diffusion timescales can 
vary significantly with changing atmospheric parameters, 
the ratios of timescales between any two elements are 
much less sensitive. 

4.2. Bulk Composition and Hydrogen 

The derived quantities (abundances, total amount of 
material found in the helium convection zone, settling 
timescales, steady state ratios and accretion rates) pre- 
sented in Table 2 reveal a number of interesting proper- 
ties about the body that polluted J0738-|-1835's photo- 
sphere. Not surprisingly, the total amount of heavy ele- 
ments is found to be mainly concentrated in the four ma- 
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Fig. 12. — Same as pre vious Figure. T he last five panels display MagE spectroscopic data. The double-peaked gas disk emission features 
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Fig. 13. — Display of our final solution ( in red ) over the lower 
resolution MMT data used in lDufour et al.l 1)20101 . see their Figure 
3). The strongest lines are identified while all the other numerous 
non-identified lines are from iron. Th e discrepancies in Fe line 
strengths noted bv lDufour et al] II2010I ) have almost all disappear 
(see text). 

jor constituent of large rocky bodies like Earth (O, Mg, 
Si and Fe) with proportions similar to those found for the 
most polluted and well studied helium- rich white dwarfs 
GD 3 62, GD 40 and HS225 3+8023 (jZuckerman et al.l 
[2007t IKlein et aLllMol [20Tll) . In total, there is some 



Wavelength (k) 

Fig. 14. — Display of our final solution (in red) over new unpub- 
lished medium resolution MMT data. Ticks and labels identify the 
strongest features present in the spectra. 

6.6 X lO^'^ g of accreted material currently present in the 
helium convection zone, the largest amount directly de- 
tected at the photosphere of a white dwarf. The object 
that tidally disrupted near J0738-I-1835 was thus at least 
as massive, and probably even more massive, than the 
dwarf planet Ceres. Assuming a typical density of ^ 2.1 
g/cm^, the alleged planet radius was at least 400 km, 
and most probably much more considering the unknown 
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Fig. 15. — Photometric measurements in GALEX, ugriz, JHK 
and Spitzer 3.6 /im and 4.5 /im bands compared to the models of 
the star and the debris disk. 
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amount of material that i) has already sunk below the 
convection zone and ii) is still contained in the surround- 
ing debris disk. 

Inspection of Table 2 also reveals a very low amount of 
hydrogen, implying that this dwarf planet was very dry, 
with less than 1% of the total mass coming from water ice 
(note that this is an upper limit since hydrogen, the light- 
est element, can only accumulate at the su rface of the 
white dwarf with cooling age) . According to lJura fc Xul 
(|201C0 . most of the internal water ice present in a 400 km 
object should survive the late stage of evolution. The low 
amount of hydrogen found at J07384-1835's photosphere 
can thus safely be interpreted as a reflection of the low 
amount of water ice initially contained in the accreted 
object. Hence, the body responsible for J0738-|-1835's 
metal pollution must have formed inside the so-called 
"snow line" where the thermodynamic conditions are 
such that rocky planetesimals initially possess little to 
no water (a similar conclusion was also recently found 
for a large sample o f nearby helium-rich white dwarfs, 
see lJura fc Xl^[20T2h . 



4.3. Oxygen Budget and Temperature Trends 

Oxygen accounts for more than half of the material 
accreted onto the surface of J0738-1-1835. Another way 
to strengthen the rocky nature of the polluting body is 
to verify that all the accreted oxygen can by accounted 
for assuming it is c ontained in rocky oxygenated minerals 
(jKlein et al.ll2010D . In other words, if the oxygen was all 
contained in mineral oxides (such as MgO, AI2O3, Si02, 
CaO, TiOa, CraOg, MnO, FeO, FcaOg, and NiO) would 
the abundance of all the heavy elements found in the 
photosphere add up in sufficient quantity to explain the 
oxygen abundanc e? Quantitatively , this corresponds to 
(see equation 3 of lKlein et al.|[20Tol) : 

2 n{0) ^ ^ p{Z) n[0) " ^ ^ 

where the heavy element Z is contained in molecule 
Zqi^z)Op(z)- Using the values presented in Table 2, we 
calculate the steady state ratios and find that the sum 
adds up only to around ~ 0.4 depending on how we di- 
vide Fe into Fe2 03 and FeO. Taken at face value, this 
indicates an oxygen overabundance, something that is 
difficult to explain given the severe constraint we have, 
from the low amount of hydrogen, on the water content 
of the accreted body. 

However, as noted above, the oxygen abundance is 
based only on 3 individual measurements which are quite 
discrepant between one another (it is not clear at this 
point whether this is due to the relatively lower resolu- 
tion and signal-to-noise ratio of the observations or dis- 
crepancies in atomic line data/broadening parameters). 
Since one of these measurements gives an abundance of 
log 0/He = -4.1 (the other two give -3.85 and -3.61), 
which is low enough to reconcile the oxygen budget, we 
could be tempted, until more precise measurements be- 
come available, to prefer the low end value in the abun- 
dance range permitted by the uncertainty. 

Alternatively, the unusually high oxygen abundance 
could be explained by the local thermodynamical condi- 
tion of the protoplanetery disk at the distance where the 
object formed. For instance, since the equilibrium com- 
position of planetesimal embr yos vary with temp erature 
and pressure at different radii ()Bond et al.ll2010a llbl). one 
expect, in a cooler environment, to observe a decrease 
in the abundances of the most refractory elements such 
as Al and C, as well as an increase in more volatile ele- 
ments like O Q. It is thus possible that most of the oxy- 
gen was actually locked up in more complex molecules 
than those assumed above (MgAl2 03, NaAlSiaOs for ex- 
ample), which would help to bring the oxygen balance 
nearer to 1. 

Interestingly, if the accreted object formed from frac- 
tionated disk material whith different thermodynamical 
conditions than where Earth formed, then it might also 
be possible to observe a signature of this in the form 
of a correlation in the abundance ratios with condensa- 
tion temperatures. Such correlations have been searched 
for in some heavily polluted and well studied helium- 

The C/O ratio is also a very im portant quantity for deter- 
mining final planetary composition {Bond et al. 2010a). Unfortu- 
nately, only an upper limit of log C/He < —3.8 is obtained from 
the absence of the strongest C transitions. 
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rich white dwa rfs such as GD 40, PG 1225-079 and 
HS 2253+8023 (IKlein et al.ll20Tol[20Tll) . 

In the case of GD 40 and PG 1225-079, the refractory 
elements are substantially enhanced relative to Mg and 
Earth values. The detailed simulations of IBond et al.l 
(|2010af ) indicate that planets forming in hotter regions 
than where Earth formed will be rich in refractories but 
low in other elements s uch as Mg, Si, O and Fe . For 
the objects analyzed by IKlein et all (|2010l 120 lit) , this 
is observed for Si, but not for Mg and Fe, suggesting 
that other processes could also be important. For in- 
stance, there could have been intense heating in the red 
giant phase, silicate vaporization, or even crust removal 
that could have affected the final bulk abundances in 
the bodies that accreted onto these stars. Consequently, 
although the planetary compositions derived in these 
stars suggest close formation to the host star, explaining 
the observed abundance pattern as a whole still remains 
quite challenging. 

With 14 heavy elements detected and measured in 
its photosphere, J0738+1835 provides us with a unique 
opportunity to empirically observe in greater details 
the abundance p attern for another large terrestrial-like 
planet. Following IKlein et al.l (|2011l . see their Figure 18), 
we plot in Figure 17 the observed abundances, normal- 
ized to Mg and bulk Earth, as a function of the conden- 
sation temperatures of each element (taking Fe, instead 
of Mg, as a reference produces a very similar plot while 
using Si shifts the points upward by about a factor of 
1.5). Our measurements suggest that there is indeed a 
correlation of the abundance ratios with condensation 
temperature, highly refractory elements being depleted 
(relative to bulk Earth) while the most volatile elements 
appear enhanced. We note that the presence of this trend 
with condensation temperature depends little on whether 
we assume that the accretion occurs in the early phase 
or at steady state (in which case a small correction, con- 
sisting of simply multiplying the results by the ratios 
of the diffusion times of the considered elements, is ap- 
plied). The pattern for the most refractory elements is 
completely the opposite to what was found for GD 40 
and PG 1225—079, clearly indicating a different forma- 
tion history and/or subsequent evolution of the body re- 
sponsible for J0738-|-1835's metal pollution. 

Rec ent simulations of th e formation of terrestrial plan- 
ets by iBond et al.l (|2010a|) suggest that objects forming 
in the inner regions have high Mg/Si, Al/Si and Ca/Si 
ratios with a steady transition toward Earth's values as 
one moves further out. This trend is a reflection of the 
condensation of refractory species in the inner most re- 
gion and of Mg-silicate further out. The low abundances 
that we find for many refractory species could thus indi- 
cate that the planet /asteroid that polluted J0738-I-1835 
formed further out in cooler regions where the Ca- and 
Al-rich inclusions (CAI), as well as other refractory-rich 
material, cannot condensate out easily. 

Although this simple interpretation is appealing to ex- 
plain our res ults, a few problern s remain. For instance, 
according to IBond et al.l (|2010al )'s simulations, while Al 
is condensing out in the inner regions of the disk, Na 
and Al are expected to condense further out together 
(species such as MgAl203 and NaAlSisOg condense out 
over a similar temperature range as Mg silicates). Hence, 
it is not clear how it could be possible to deplete Al and 
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Si but not Na (a similar argument can be applied to Mn 
as well). 

It must be noted, however, that the error bars on some 
of these elements are quite large, severely affecting our 
level of confidence on the existence of the trend for the 
low range of condensation temperature (Tc < 1200 K). 
In particular, Na and Mn abundances rely on very few 
measurements and given the signal-to-noise ratio of our 
observations in the region where these lines are found 
(and the corresponding error bars), it is quite possible 
that the real abundances of these elements could lie on 
the low end range of our measurements, which would be 
more compatible with Earth-like values. Furthermore, 
it must be noted that the estimation of the composition 
of b ulk Earth volatile mater ial is quite a difficult task 
(see lAllegre et al.lll995l [200ll ) and that uncertainties on 
the Na and Mn Earth values used for the normalization 
in Figure 17 are to be considered as well. For exam- 
ple, taking N a and Mn values from lAllegre et al.l (jl995l ) 
over those of lAUegre et al.l (|200lD would lead to a bet- 
ter agreement between the observed abundances of these 
elements and Earth values. 

Since, as we discussed above, the oxygen abundance 
is also quite uncertain (the lower abundance of oxygen 
would also be more compatible with the oxygen balance 
argument in certain conditions), another coherent pic- 
ture could be that the observed abundances are in fact 
close to Earth values for most elements, with only a few 
exceptions like Si, Ti and Al, which can confidently be 
considered depleted relative to bulk Earth. 

Such a pattern, if confirmed with better accuracy, 
could be interpreted very differently. For example, it is 
possible that differentiation took place in the body that 
accreted at the surface of J0738+1835 (in fact, with a 
minimal radius of at least 400 km, diffe rentiatio n most 
probably occurred). As discussed in Klei n et "all ([2010), 
elements such as Si, Ga, Ti, Al tend to concentrate in 
the crust of a differentiated body. Given that it is these 
very same elements that are found to be depleted in our 
object, the abundance pattern found in J0738-f 1835's 
photosphere could very well be due to some sort of crust 
removal. A similar process may explain the strange abun- 
dance pattern found in GALEX J1931-^0117 ((Melis et all 
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120111) . where the planet's outer layers is believed to have 
been stripped away by some wind interaction in the AGB 
phase. Another possibility would be that a large impact 
striped away a large portion of the crust, affecting the 
final composition that we observe at th e white dwarf 
photosphere (see iZuckerman et all 120 111 where a simi- 
lar scenario has been proposed to explained the peculiar 
abundances found in NLTT 43806). 

It is unfortunate that a detailed analysis such as this 
one, even using data from some of the best observational 
facilities in the world, cannot discriminate furthermore 
between the different possible scenarios. Although our 
analysis clearly indicates that the most refractory ele- 
ments are depleted relative to Mg and Earth, the case 
for a trend with condensation temperature must be con- 
sidered more speculative at this point due to the uncer- 
tainties in our measurements. . 

Future investigation aiming at refining, among other 
things, the oxygen, sodium and manganese abundances, 
should help to rule out some of the various scenarios de- 
scribed above. Meanwhile, until better spectroscopic ob- 
servations become available, we prefer to remain cautious 
and not overinterpret the data. Regardless, our study 
seems to indicate that that we are on the verge of enter- 
ing an exciting new era of high precision measurements of 
terrestrial-like planet compositions from polluted white 
dwarfs and that appreciable insights on planetary forma- 
tion should be accessible in the near future.. 

5. CONCLUSION 

We presented a detailed analysis, based on high res- 
olution optical spectroscopic observations, of the DBZ 
SDSS J073842. 56-^183509.06, the most metal polluted 
white dwarf currently known. We also obtained IRAC 
Spitzer infrared photometry at 3.6 fim and 4.5 /im, 
which, combined with earlier JHK photometric mea- 
surements, allows us to better constrain the properties 
(inner and outer temperature, inclination) of the bright 
debris disk surrounding that white dwarf. 

We have measured the abundances of 14 different 
heavy elements and determined that the total mass of 
material currently present in the helium convection zone 
is ~ 6.6 X lO^'^ g, indicating that the object that was 
tidally disrupted was at least as large as the dwarf planet 
Ceres. We also find, from the low amount of hydrogen 
present in the photosphere, that the total amount of wa- 
ter ice was less the 1% of the mass of this object, sug- 
gesting that it formed inside the so-called snow line. 

Our most interesting finding, summarized in Figure 17, 
is that the highly refractory species are depleted, relative 
to Mg and Earth, while more volatile elements appear to 
be enhanced, with a possible correlation with the con- 
densation temperature. Such a trend could be the signa- 
ture that the dwarf planet that accreted on J0738-(-1835 
formed in a lower temperature environment than Earth, 
where the thermodynamical conditions are such that the 



planetesimal embryos' chemical composition tends to be 
depleted of more refractory elements. However, the cor- 
relation must be considered uncertain given the large un- 
certainties associated with the O, Na and Mn measure- 
ments. Nevertheless, it is clear that the Si/Mg, Ti/Mg 
and Al/Mg ratios are lower than bulk Earth values, indi- 
cating that the physical processes that affected the final 
composition of the dwarf planet were different than those 
that occurred for Earth and other polluted white dwarfs 
such as GD 40, GD 362 and PG 1225-079. 

Another way to interpret our findings, if we momentar- 
ily ignore the more uncertain O, Na and Mn abundances, 
is that differentiation occurred in the large body that 
produced the pollution observed on J0738-I-1835. With 
a radius of at least 400 km, a crust containing prefer- 
entially Si, Ti, Al and Ca could have formed and been 
removed, potentially explaining part of the abundance 
pattern that we find. 

Improvements, allowing us to discriminate between 
these various scenarios, albeit difficult, should be ob- 
tainable from longer observing runs on large telescopes 
(the combination of high resolution with high signal-to- 
noise ratio is the key to obtain precise results). This 
work, as well a many other recen t studies (for ex- 
ample. IZuckerman et ahl 120071 l2010t iMelis et al.l 120111: 
IKlein et al.l I2010L 120111) demonstrate, once more, the 
extraordinary potential of metal contaminated white 
dwarfs to study rocky extrasolar planet /asteroid bulk 
compositions. These objects are unique probes that can 
serve as perfect testbeds for planetary formation and evo- 
lution theories, providing information that cannot be ob- 
tained from any other way. As our understanding of 
the physical processes at play increases, as well as our 
confidence in the interpretation of the chemical patterns 
found in these stars, reliable bulk composition of rocky 
extrasolar objects, conceivably with a level of precision 
comparable, if not better, to what is possible for simi- 
lar objects in the solar system, should soon be provided 
from metal polluted white dwarf studies. 
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